Bose-Einstein condensation of 86 Sr 
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We report on the attainment of Bose-Einstein condensation of 86 Sr. This isotope has a scattering 
length of about +800 ao and thus suffers from fast three-body losses. To avoid detrimental atom 
loss, evaporative cooling is performed at low densities around 3 x 10 12 cm -3 in a large volume optical 
dipole trap. We obtain almost pure condensates of 5 x 10 3 atoms. 

PACS numbers: 37.10.De, 67.85.-d, 67.85.Hj 



Quantum degenerate gases of atoms with two valence 
electrons are an exciting field of research. The elec- 
tronic structure of these atoms is the basis for appli- 
cations like ultraprecise optical clocks [H-Q- Quantum 
degenerate samples open new possibilities such as the 
study of unique quantum many-body phenomena [n- 
[llj and novel schemes of quantum computation [l2l-[l5j. 
They are also an ideal starting poin t for the creation of 
molecules of two-electron atoms \Ya a nd molecules of 
two-electron atoms with alkali atoms jl&[l9l , which have 
applications in precision measurement [20j or quantum 
simulations with spin-dependent anisotropic long range 
interactions [21[. The first two-electron atom cooled to 
quantum degeneracy was ytterbium [22l-[24j , followed by 
calcium [25| and strontium last year. Quantum degener- 
ate samples of bosonic 84 Sr and 88 Sr [26l428j, as well as 
of fermionic 87 Sr [29|, [30| have been created, but Bose- 
Einstein condensation (BEC) of 86 Sr has so far been elu- 
sive [ail. 

The success of evaporative cooling with the goal to 
reach quantum degeneracy depends largely on the ratio 
of elastic to inelastic collisions of the atomic species used 
(32j . The four isotopes of strontium differ significantly in 
this respect. Despite its low natural abundance of 0.56%, 
84 Sr was the first to be Bose condensed as it offers ideal 
scattering properties [26|, [27j . 88 Sr, which has a scat- 
tering length close to zero, and spin-polarized fermionic 
87 Sr barely collide at ultralow temperatures. They could 
only be cooled to quantum degeneracy using mixtures 
with other isotopes or spin states [28-130] . 86 Sr poses the 
opposite challenge: its large scattering length of +800 ao 
[33( | leads to a large three-body loss rate constant [3l| . 
Previously reported conditions of 86 Sr in a dipole trap 
led to an elastic to inelastic collision ratio that was insuf- 
ficient for evaporative cooling [31]. Magnetic tuning of 
the scattering properties, as was an essential ingredient 
to circumvent similar difficulties on the way to BEC of 
cesium (34|, is not available in alkaline-earth elements. 
Another approach to increase the ratio of elastic two- 
body collisions and inelastic three-body collisions is to 
lower the density. This approach was important for BEC 
of 40 Ca, which is also an isotope with large scattering 
length [25j, as well as for BEC of cesium. 

In this Rapid Communication, we show that it is 
indeed possible to create 86 Sr BECs by performing 



evaporation at comparatively low densities around 3 x 
10 12 cm -3 . To obtain sizeable atom numbers at these 
low densities, we implement a large volume optical dipole 
trap. 

Our experimental procedure initially follows the one 
used for Bose condensation of 84 Sr [26j. Zeeman slowed 
atoms are captured and cooled by a "blue" magneto- 
optical trap (MOT) operating on the 5s 2 1 5'o — 5sbp 1 P\ 
transition at a wavelength of 461 nm. A weak leak of 
the excited state of the cooling cycle continuously pop- 
ulates the 5s5p 3 P2 metastable state. Weak-field seek- 
ing atoms in this state can be trapped in a magnetic 
trap formed by the quadrupole magnetic field used for 
the MOT. Metastable state atoms are accumulated in 
the magnetic trap for Is. Further cooling and den- 
sity increase is achieved by operating a "red MOT" on 
the 7.4 kHz linewidth 1 So — 3 Pi intercombination line 
at 689 nm. To increase the capture velocity of the red 
MOT, we frequency modulate the light, producing side- 
bands, which cover a detuning range from a few ten kHz 
to a few MHz to the red of the transition. To load the 
red MOT, the metastable-state atoms in the magnetic 
trap are optically pumped to the 1 So ground state us- 
ing the 5s5p 3 P2 — 5s5d 3 D2 transition at 497 nm. After 
loading, the MOT is compressed by reducing the MOT 
beam intensity and ramping off the frequency modula- 
tion, resulting in a colder and denser sample. At this 
point the MOT contains 2 x 10 T atoms at a temperature 
of 1 fiK. 

The following evaporative cooling stage significantly 
differs from our previous experiments on 84 Sr. The 
atoms are transferred into a crossed-beam optical dipole 
trap, which is much larger than the one used before. The 
trapping geometry consists of a horizontal and a nearly 
vertical beam, derived from a broadband ytterbium fiber 
laser operating at 1075 nm. The horizontal beam has an 
elliptic beam shape with a horizontal waist of 300 /im and 
a vertical waist of 33 /im. The vertical beam is circular 
with a waist of 290 /im and is used to provide additional 
confinement along the weak direction of the horizontal 
beam towards the end of evaporation. Initially, only the 
horizontal dipole trap beam is used and set to a power 
of 2 W. The resulting trap is oblate with horizontal trap 
oscillation frequencies of 3 Hz and 30 Hz, a vertical trap 
oscillation frequency of 260 Hz and a potential depth of 
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3.7 /iK, taking into account gravitational sagging. This 
oblate trap geometry combines the requirement of strong 
enough confinement against gravity in the vertical direc- 
tion with the requirement of a large trap volume. 

For optimum loading of the dipole trap, we adjust the 
intensity and detuning of the red MOT beams before 
switching them off. After 1 s of plain evaporation, 2.5 x 
10 6 atoms at a temperature of 500 nK remain in the trap. 
The peak density is 3 x 10 12 cm -3 and the peak phase- 
space density is 0.05. Even at the temperature of 500 nK, 
unitary limitation leads to a reduction of the thermally 
averaged elastic scattering cross section by a factor of 
two compared to the zero temperature value. The elastic 
collision rate is 380 s _1 . 

After the plain evaporation stage, forced evaporation 
is performed over 4.8 s. The power of the horizontal 
dipole trap beam is reduced nearly exponentially to 
740 mW, with a longer time constant during the last 
1.5 s. During the first 2.3 s of forced evaporation, the 
vertical dipole trap beam is increased to a power of 
0.5 W, at which it stays for the rest of the evaporation 
sequence. The vertical beam increases the confinement 
along the weak direction to 5 Hz and does not provide 
any confinement along the vertical direction. The po- 
tential depth of the vertical beam in the radial direc- 
tions is 200 nK. During evaporation atoms escape the 
trap mainly downwards by leaving the trap over the po- 
tential barrier formed by the horizontal beam and grav- 
ity. It is beneficial for evaporation that the vertical trap 
oscillation frequency remains higher than the elastic col- 
lision rate, so that high energy atoms produced by col- 
lisions can quickly move out of the trap. During evapo- 
ration the peak density slightly drops to 10 12 cm -3 and 
the elastic collision rate drops to 70s -1 before conden- 
sation. At this point the temperature is ~30 nK and the 
thermally averaged scattering cross section is very close 
to the zero- temperature value. 

The phase transition from a thermal cloud to a BEC 
becomes evident in the appearance of a bimodal dis- 
tribution, as clearly visible in time-of-flight absorption 
images and in the corresponding linear density profiles 
shown in Fig. [TJ At higher temperature the distribution 
is thermal, exhibiting a Gaussian shape. Cooling below 
the critical temperature T c leads to the appearance of an 
additional, narrower and denser, elliptically shaped com- 
ponent: the BEC. The phase transition occurs after 4 s 
of evaporative cooling, when the power of the horizontal 
beam is 820 mW. At this point, the horizontal trap os- 
cillation frequencies are 5 Hz and 18 Hz, the vertical trap 
oscillation frequency is / vert = 115 Hz and the trap depth 
is about 110 nK, taking into account gravitational sag- 
ging. Note that despite the large scattering length, we 
are not in the collisionally hydrodynamic regime since 
the elastic scattering rate is only one tenth of 27r/ ver t 

m. 

In order to analyze the phase-transition precisely, we 
fit bimodal distributions to the time-of-flight absorption 
images. The bimodal distributions consist of a Gaus- 
sian for the thermal part of the cloud and an inverted 
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Figure 1: (Color online) Absorption images and integrated 
density profiles showing the BEC phase transition for differ- 
ent times t ev of the evaporative cooling ramp. The images 
are taken in the horizontal direction at 45° to the horizon- 
tal trap axes, 25 ms after release from the trap. The field 
of view is 400 /mi by 400 /xm. The solid line represents a fit 
with a bimodal distribution, while the dashed line shows the 
Gaussian-shaped thermal part, from which the given temper- 
atures are derived. 



paraboloid integrated along one direction for the BEC. 
We extract the amount of condensed atoms No , the total 
amount of atoms N and the temperature T of the sam- 
ple from the fits [36j . The condensate fraction Nq/N as 
a function of temperature, scaled by the transition tem- 
perature of a non-interacting gas in the thermodynamic 
limit is shown in Fig. [2j From this data we extract a 
critical temperature of T c = 29(2) nK= 0.90(5) T c °. The 
atom number at the onset of BEC is 3.5 x 10 4 . Inter- 
action effects reduce the critical temperature T c by 14% 
compared to Tj, whereas finite-size effects add only a 
minor reduction of 1%. These effects give a theoreti cally 
expected transition temperature of T c = 0.85 T® [37[, 
which is within the uncertainty of our data. 

After 4.8 s of evaporation, a nearly pure BEC of 5 x 10 3 
atoms is produced in a trap with oscillation frequencies 
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Figure 2: (Color online) The condensate fraction No/N as a 
function of the temperature scaled by the critical temperature 
of an ideal gas in the thermodynamic limit T<?. The line 
shows the behavior expected for an ideal gas. The arrow 
indicates the expected transition temperature T c = 0.85 
when interactions and minor contributions from finite size 
effects are taken into account. 
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Figure 3: (Color online) Inversion of aspect ratio during the 
expansion of a pure BEC. The images (fields of view 300 /mi 
x 300 /im) are taken along the horizontal direction and show 
the BEC in the trap (left) and after 25 ms free expansion 
(right). 

of 5 Hz and 16 Hz in the horizontal plane and 70 Hz in 
the vertical direction. The peak-density of the BEC 
is 3 x 10 12 cm -3 and the chemical potential is 7nK. 
The BEC has an oblate shape with calculated in-situ 
Thomas- Fermi radii of 40 jiva and 12 /xm in the horizon- 
tal plane and 3 /im in the vertical direction. After release 
from the trap, the BEC undergoes a mean-field driven 
expansion, with the strongest expansion in the vertical 
direction, resulting in an inversion of ellipticity as shown 
on the absorption images in Fig. [3l 

We extract an upper bound for the three-body loss 
rate constant from the lifetime measurement of a nearly 



pure BEC. The initial atom loss rate per atom from the 
BEC is N^dNo/dt = 3(1) s" 1 . This value is determined 
from the initial slope of an exponential fit to 500 ms of 
decay data. Assuming only three-body loss, the loss 
of atoms is described by dNo/dt = — (1/6) if 3 J n^dV, 
where no is the density distribution of the BEC and 
the factor 1/6 takes into account the difference of three- 
body correlations of a thermal gas and a BEC. The re- 
sulting upper bound for the three-body loss rate con- 
stant is K3 = 6(3) x 10 _24 cm 6 /s. This value is about 
six times higher than the value measured by Ferrari et 
al [3l|| . Moreover it is an order of magnitude higher 
than the maximally expected loss rate constant, which is 
K3 = 210ha 4 /m = 5 x 10 -25 cm 6 /s, assuming that three 
atoms are lost during each three-body loss event [38|, [39| . 
Here a is the scattering length and m is the mass of 86 Sr. 
Several explanations are possible for the unusually large 
measured value. More than three atoms might be lost 
per three-body loss event resulting from secondary col- 
lisions, possibly augmented by an enhanced atom-dimer 
scattering cross section [Io|, El| . Evaporation of atoms 
heated by technical noise of the dipole trap is an example 
for an explanation other than three-body loss. Further 
studies are needed to resolve this issue. 

In conclusion, we have produced a Bose-Einstein con- 
densate of 86 Sr containing 5 x 10 3 atoms. With this 
achievement all stable isotopes of strontium have been 
cooled to quantum degeneracy. The large scattering 
length of 86 Sr leads to a very large three-body loss rate 
coefficient, which poses a challenge for evaporative cool- 
ing. We have shown that performing evaporation at low 
density in a large volume trap is a possible way to over- 
come this problem. The BEC of 86 Sr enriches the pos- 
sibilities opened up by strontium quantum gases. The 
large scattering length of 86 Sr originates from a weakly 
bound state in the molecular potential. This state will 
influence the properties of optical Feshbach resonances 
H2,IIj|. The BEC of 86 Sr increases the options for quan- 
tum degenerate mixtures of different Sr isotopes or of 
mixtures of Sr with other elements. This wider choice 
can be important, for example in the search for mag- 
netic Feshbach resonance in Sr-alkali mixtures suitable 
for molecule creation [l8| . 
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